Coreless vortex formation in a spinor Bose-Einstein condensate 
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Coreless vortices were phase-imprinted in a spinor Bose-Einstein condensate. The three- 
component order parameter of F = 1 sodium condensates held in a loffe-Pritchard magnetic trap 
was manipulated by adiabatically reducing the magnetic bias field along the trap axis to zero. This 
distributed the condensate population across its three spin states and created a spin texture. Each 
spin state acquired a different phase winding which caused the spin components to separate radially. 
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Spin textures play a central role in describing the 
physics of elementary particles jl), liquid '^He-A [^,^ 
the quantum Hall effect |5|], and gaseous Bose-Einstein 
condensates [||, 0, ^, ||, ^|. Topological defects vary be- 
tween superfluid systems described by scalar and vector 
order parameters. In spin-less or spin-polarized conden- 
sates, line defects such as vortices have cores where the 
density of condensed particles is necessarily zero to keep 
the order parameter sing le- valued |n|, |l|, |l|. However, 
in condensates with an internal, spin degree-of-freedom, 
coreless vortices exist as spin textures [10. Such struc- 
tures are referred to as skyrmions (Anderson- Toulouse 
vortices |3|) or merons (half-skyrmions, Mermin-Ho vor- 
tices [^) depending on the boundary conditions of the 
system. 

In this Letter, we study spin textures in a Bose- 
Einstein condensate (BEC). Coreless vortices were cre- 
ated in F = 1 spinor condensates held in a loffe-Pritchard 
magnetic trap by adiabatically reducing the magnetic 
bias field along the trap axis to zero. This continuously 
transformed the initially spin-polarized condensate into 
a coherent superposition of three spin states, each with 
a different phase winding. The resulting angular mo- 
mentum per particle varied between spin states and the 
condensate evolved such that states with more angular 
momentum per particle circulated around states with less 
angular momentum per particle. Thus, the condensate 
had a net axial magnetization that varied with radial po- 
sition. Previous work on vortices in a two-component 
system used laser, microwave, and radio frequency fields 
to spatially and temporally control the interconversion 
between components jl^. However, without these ap- 
plied fields the two-components evolved independently as 
distinguishable fluids. In our work, the spin states can 
freely interconvert at all points in space and time such 
that the spin texture would continually heal itself even 
in the presence of state-dependent losses. 

In cylindrical coordinates, the spin— _F conden- 
sate wavefunction can be written as |^'(r, z)) = 
\Jn(r, (f), z)|C(r, cj), z)), where n is the atomic number 
density and the 2F + 1 component spinor |C) = 
J2m,=-FCmjF,mz), |(C|C)P = 1 describes a spin tex- 



ture. A loffe-Pritchard magnetic trap consists of an 
axial bias field (with curvature) and a two-dimensional 
quadrupole field in the orthogonal plane |l6| : 

B{r, 0, z) = BJ + B'r{cos{2(j)) f - sm{2(j)) $), (1) 

where B' is the radial magnetic field gradient and 
quadratic terms have been neglected. For a condensate 
of radial extent R confined in a loffe-Pritchard magnetic 
trap with Bz B'R > 0, |C} — \F,mz = mp), where 
niz and mp are the projection of the atomic spin along 
the z-axis and local magnetic field direction, respectively. 
Adiabatically ramping Bz from Bz ^ B'R > to zero 
rotates the atomic spin about the position-dependent 
axis n{(l)) — sincj) x + coscj) y, and drives the tran- 
sition |F, = mp) I]ml=-|mp| Cm. exp(«(m2 - 
■mF)(f) \F,mz) [0 Thus, the condensate popula- 

tion is distributed across 2|r7ii?| -I- 1 spin states with each 
acquiring a different topological phase factor and angu- 
lar momentum per particle due to the variation of Berry's 
phase with magnetic quantum number |l9| . 

The condensate remains in the state \F,mp) with re- 
spect to the local magnetic field provided the local Zee- 
man energy, ~ gp^B{B^ + (i3'r)^)^/^, dominates the 
local kinetic energy associated with the spin texture, 
~ 'h? /mr"^ , where gp is the Lande g-factor, /is is the Bohr 
magneton, and m is the atomic mass. For Bz = 0, atomic 
spins aligned with the quadrupole magnetic field produce 
the planar spin texture in Fig. 0(a). However, the infi- 
nite kinetic energy associated with the wavefunction sin- 
gularity at r = creates a non-planar spin texture over a 
disc of radius {h^ /mgp^BB'Y^^ , with higher angular 
momentum spin states residing outside those with lower 
angular momentum. 

Bose-Einstein condensates containing over 10^ ^^Na 
atoms were created in the 1) state in a magnetic 
trap, captured in the focus of an optical tweezers laser 
beam, and transferred into an auxiliary "science" cham- 
ber as described in Ref. j2^. In the science cham- 
ber, the condensate was loaded into a microfabricated 
loffe-Pritchard magnetic trap formed by a Z-shaped 
50 /im X 10 /im electroplated copper wire carrying cur- 
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(a) 








Along the wiretrap axis, the magnetic field was 



FIG. 1: Planar spin textures. Spins aligned with two- 
dimensional (a) quadrupole and (b) hexapole magnetic fields 
produce textures with winding numbers —1 and —2, respec- 
tively. In general, spins aligned with a two-dimensional ax- 
isymmetric multipole magnetic field, B = Ar'^^^ {cos{q(f)) f — 
sm{q(f)) (p), q > 0, produce a texture with winding number 
— (g — 1). Counter-clockwise traversal of the dashed contours 
in (a) and (b) leads to clockwise (negative) spin rotation, 
with the winding number defined as the integer number of 
revolutions made by the spin vector while circumnavigating 
the singularity at the origin ^]. 



rent / and an external magnetic bias field, B±, as de- 
scribed in Ref. | |2^ . Typical wiretrap parameters were 
/ = 720 mA, Bj_ = 5.3 G, and = 1.3 G, resulting in a 
radial magnetic field gradient B' = 180 G/cm. This pro- 
duced axial and radial trap frequencies lUz — 2tt x A Hz 
and uj± — 2n X 250 Hz, respectively. Condensates in the 
wiretrap had > 2 x 10^ atoms, a Thomas-Fermi radius of 
« 5 /im, and a lifetime > 25 s 
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Coreless vortices imprinted onto the condensate wave- 
function by adiabatically ramping B^ ^ are shown in 
Fig. ^ To observe the nature of the spin texture, an axial 
bias field was switched on non-adiabatically along either 
the negative [Fig. ||(a)-(d)] or positive [Fig. ||(e)-(h)] z- 
axis. Switching the axial bias field on suddenly "froze" 
the atomic spins and effectively "projected" the conden- 
sate wavefunction onto a basis quantized with respect to 
the local (axial) magnetic field. This allowed the spin 
states to be separated by a magnetic field gradient ap- 
plied during ballistic expansion. Switching the direction 
of the axial projection field exchanged the roles of the 
|1,— 1) and |1,+1) states. Figures ||(a) and |(e) show 
the coreless nature of the vortices, while Figs. |(b) and 
^(f) show the concentric cylinder structure resulting from 
the competition between the atomic Zeeman energy and 
the kinetic energy of the rotating spin states. We assume 
that the two-dimensional (cj^ <^ ballistic expansion 
process simply magnifies the condensate wavefunction, as 
it does in the expansion of a single-component conden- 
sate with vortices [|4[ p5| . 

To imprint the coreless vortices, the axial bias field was 
ramped linearly from Bz = 1.3 G to i?z « in 10 ms. 



B{r = 0,cj),z) = {Bz 



(2) 



where quadratic terms neglected in Eq. (|^) are included. 
The presence of axial magnetic field curvature, B" k, 
5 G/cm^, implies that the spin texture had a slight axial 
dependence. Ramping Bz ^ Q quickly (10 ms) com- 
pared to the axial trap period (250 ms) compressed the 
condensate radially but left the axial extent of the con- 
densate unchanged. Thus, the axial magnetic field varia- 
tion, Ai3z, was less than the radial magnetic field varia- 
tion, ABr, and the axial dependence of the spin texture 
was minimal. With Bz = 1.3 G [Bz « G), the con- 
densate chemical potential was /i ~ (/iB/2) x 3 mG (/i ~ 
(/ii3/2) X 27 mG) yielding AB^ « 3 mG (AB^ « 27 mG). 
The images shown in Fig. |^ integrated the atomic num- 
ber density along the z-axis and therefore averaged over 
the minor axial variation to the spin texture. 

To project the condensate wavefunction onto a basis 
quantized with respect to the local magnetic field, an ax- 
ial bias field was switched on at a rate ofi?z = 2x10^ G/s 
along either the negative or positive z-axis. For < Br < 
ABr, the Landau-Zener non-adiabatic transition proba- 
bility, exp{—'KfiBBr/hBz) > 0.9, was sufficiently close to 
unity that the atomic spins remained "frozen" during the 
sudden application of the axial bias field and the spin 
texture could be accurately diagnosed. The magnetic 
trap was switched off 100 /iS after applying the axial bias 
field allowing the atoms to expand ballistically. While 
the total condensate density monotonically decreased as 
a function of radial position [Figs. ||(a) and ||(e)], the 
density of each spin component peaked at a different ra- 
dius signifying a variation in the angular momentum per 
particle between spin states [Figs. ||(b) and ||(f)]. 

With the application of a projection field along the 
positive z-axis, additional rings of atoms appeared in the 
1,-1) and 11,0) states [Fig. ||(f)]. Possible technical 
origins include the non-instantaneous turn on of the pro- 
jection field and the ballistic expansion process. Also, 
ramping — > in 10 ms was non-adiabatic for atoms 
near r — causing an atom loss of w 50% as observed 
in previous work [Q. If these atoms had not left the 
condensate before the projection field was applied, they 
may have contributed to the images displayed in Fig. ^(f). 
The additional rings of atoms may also correspond to a 
low energy, radial spin- wave excitation ||, ^. However, 
we could not identify any asymmetry between applying 
the projection field along the positive versus negative z- 
axis that would account for the presence of the extra rings 
in Fig. |(f) but not in Fig. |(b). 

Engineering topological states in a Bose-Einstein con- 
densate has received much theoretical attention []l^ |2^, 
p7| , p8| , p9| |30[ |. The evolution of a condensate confined 
in a loffe-Pritchard magnetic trap while ramping Bz 
is described by a position-dependent spin rotation about 




FIG. 2: Coreless vortex formation in a spinor Bose-Einstein condensate. Coreless vortices were imprinted by ramping 
and diagnosed by suddenly switching (a)-(d) Bz ^ and (e)-(h) Bz ^ 0. Axial absorption images display the optical density 
of condensates after 20 ms of ballistic expansion (a,e) without and (b,f) with a magnetic field gradient applied to separate the 
different spin states, (c) and (g) Azimuthally averaged optical density vs radial position for spin components shown in (b) and 
(f), respectively. The radial separation of the spin states resulted from their relative phase windings and is a clear signature of the 
skyrmion/meron wavefunction [Eq. (^]. (d,h) Axial magnetization per particle, M/N = {fiB/'2) x {hmp=+i — nmp=-i) /fitotai, 
vs radial position. The absorption imaging light was resonant with the _F = 2 ^ F' = 3 transition. The atoms were optically 
pumped into the F — 2 hyperfine level with a pulse resonant with the F = 1 F' — 2 transition. This provided equal imaging 
sensitivity to each spin state. The field of view in (a), (b), (e), and (f) is 1.0 mm x 3.0 mm. 



the n{(j)) axis through an angle (3{r) [|T8| 

|C(r,^,z))=e-^(^/'')-"(«'5M|Co), (3) 

where T is the spin operator and |Co) — \F,mz ~ mp) is 
a polarized spinor. For |Co) = 1); Eq. (||) gives the 
condensate spinor in the laboratory frame as 

cos2(/3(r)/2) |1,-1) + 
|C(r,0,z))= ^sin(/?(r)) e^^ |1,0)+ (4) 
sin2(/3(r)/2) e*^^ |1,+1). 

/3(0) — and skyrmions (merons) are described by the 
boundary condition /3(oo) = tt (/3(oo) = tt/2). 

The radial dependence of /3(r) is determined by requir- 



ing the Gross-Pitaevskii energy functional, 

(5) 

be stationary with respect to variations in (3. In writing 
Eq. we have made the Thomas-Fermi approximation, 
(VCIVCH El=-F Vd. ■ VC™., V = -gpf^BF ■ B/h, 
F — Co = 'inh'^a/m, and C2 = inh'^Aa/m. Here 

a = (2ao -l- Q2)/3 and Aa = (02 — ao)/3 characterize two- 
body interactions, where oq and a2 are scattering lengths 
for collisions with total angular momentum _F = and 
F = 2, respectively [^. For 02 > the atomic inter- 
actions are anti- ferromagnetic (polar), as in ^^Na [^ , 
while for 02 < they are ferromagnetic, as in ^^Rb [ p2[ . 

Using Eq. (||) as a trial spinor, we find that (3{r) sat- 
isfies the meron boundary conditions and varies from 
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/3(0) = 0, due to the kinetic energy of the spin tex- 
ture, to /3(oo) = 7r/2, due to the atomic Zeeman energy, 
over a characteristic length scale given by the larger of 
(fi^ /mgp^BB'Y^^ and \Bz\/ B' , in reasonable agreement 
with the results presented in Fig. ^j. We observed that 
the boundary condition /3(0) — was maintained regard- 
less of the sign of Bz, i.e. atomic spins along the trap axis 
always remained in the initial state. This emphasizes the 
non-trivial topology of the spin texture. If the atomic 
spins had simply followed the local magnetic field, (3{r) 
would satisfy tan/3(r) = B'r/Bz- Thus, scanning B^ 
from slightly positive to slightly negative would instan- 
taneous change the boundary condition at the origin from 
/?(0) = to (3{0) = TT and flip the atomic spins along the 
trap axis. 

The trial spinor in Eq. does not have the most 
general form since it was derived by rotating a polar- 
ized spinor [Eq. (^] and is inherently ferromagnetic, 
|i^| — h. Accordingly, the spin-dependent interaction 
term in Eq. does not vary with (3 and therefore does 
not contribute to the determination of /3(r). This restric- 
tion is not severe since the Zeeman energy dominates the 
spin-dependent interaction energy in our experiment. 

It should be possible to overlap the condensate with an 
optical dipole trap so that, after ramping Bz 0, relax- 
ing the radial magnetic field gradient would allow for the 
spin-dependent interactions to determine the evolution of 
the condensate. It is predicted that skyrmions/merons in 
condensates with anti- ferromagnetic (ferromagnetic) in- 
teractions are unstable (stable) [^ |l^ |3[ |^ |35[ |3^, 
|37| , At present, excitations created during the field 
ramping process have prevented us from characterizing 
the stability of coreless vortices in an anti-ferromagnetic 
^^Na condensate. However, we were able to imprint the 
spin texture in the presence of the optical dipole trap, as 
well as reproduce the results of Ref. [^ by fully inverting 
the axial bias field to produce vortices with a An phase 
winding in an optical dipole trap. At zero magnetic field, 
it may then be possible to observe multiply charged vor- 
tices in a spinor condensate "unwind" themselves as pre- 
dicted in Ref. [^ . The unwinding process is precisely the 
position-dependent spin rotation demonstrated in previ- 
ous work with the spin texture investigated here as 
an intermediate state pTf . 

In conclusion, we have demonstrated a robust tech- 
nique for creating coreless vortices in a Bose-Einstein 
condensate. Our technique can be extended to generate 
spin textures with arbitrary winding number and variable 
angular momentum per particle by using higher-order, 
axisymmetric multipole magnetic fields and condensates 
with different spin. This work opens up the opportu- 
nity to study the stability of topological defects in spinor 
Bose-Einstein condensates. 
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and Lucile Packard Foundation. 
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